Hydrodynamical evolution of dissipative QGP fluid 
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■ Abstract. In the framework of the Mueller-Israel-Stewart theory of dissipative fluid dynamics, 

we have studied the space-time evolution of a QGP fluid. For simplicity, we have considered 
shear viscosity only and neglected bulk viscosity and heat conductivity. Shear viscosity opposes 
the expansion and cooling of the fluid. As a result, the lifetime of the fluid is extended. We also 
find that the parton pr distribution is considerably flattened. 
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1. Introduction 

One of the primary goals of relativistic heavy ion collisions is the creation and detection of 
the (lattice QCD pQ) predicted deconfined phase of quarks and gluons. Experiments at RHIC, 
BNL, gave strong indications of QGP formation in central Au+Au collisions 2 J. Hydrodynamics 
provides a simple, intuitive dynamic description of relativistic heavy ion collisions. Most of the 
experimental data on soft hadron production at RHIC are well described by ideal hydrodynamics. 
However, some problems remain. For example, one of the most important signals of collectivity, 
the elliptic flow, saturates beyond px ~ 1.5—2 GeV/c while ideal hydrodynamics predicts a 
continuous increase with px [3] • 

Ideal hydrodynamic may not give an accurate description of the early stage of the collision. 
In the early stages of collision, particle momenta are predominantly in the beam direction, while 
ideal hydrodynamic assumes a locally isotropic distribution. The anisotropy in the momentum 
distribution will give rise to shear viscous stress. Other dissipative effects, e.g. bulk viscosity and 
heat conduction can also affect the hydrodynamic evolution of QGP. Unlike ideal flow, non-ideal 
flow generates entropy. The space-time evolution is also changed. 

Relativistic theories for dissipative fluid dynamics were formulated long ago [IJ 03 El • in early 
formulations, called 1st order theories HH3J) the entropy 4-current contained terms of first order 
in the dissipative fluxes; these had the undesirable feature of violating causality. This is corrected 
in so-called 2nd order theories, due to Grad, Mueller and Israel and Stewart 0. In 2nd order 
theories, the entropy 4-current contains terms of 2nd order in the dissipative fluxes. The space of 
thermodynamic variables is extended to include the dissipative fluxes, and relaxation equations 
for these dissipative fluxes are obtained from the positivity condition for entropy production 
> 0. 

Up to now, only a few authors have considered the effects of dissipation in relativistic 
heavy ion collisions El El EDS EEU El EES! • Most of these studies considered one dimensional 
systems using 1st order theories. Only recently 2nd order theories have been implemented 
numerically ^QIEIEII- In the present paper, using the causal 2nd order theory of dissipative 
fluid dynamics, we study the effect of shear viscosity on the space-time evolution of a QGP 



fluid. The paper is organized as follows: In Section II, we present the equations for 2nd order 
dissipative hydrodynamics. Our choice of equation of state, transport coefficients and initial 
conditions is described in Section III. Results for the space-time evolution of the QGP fluid, 
with and without dissipation, are presented Section IV. Conclusions are drawn in Section VI. 



2. Equations for causal dissipative hydrodynamics 

We consider a QGP fluid in the central rapidity region, with zero net baryon density and chemical 
potential, ns = 0, fis = 0. We neglect the effects of heat conduction = 0) and bulk viscosity 
(massless particles) and account only for shear viscosity. We work in the Landau-Lifshitz energy 
frame. 

The energy- momentum tensor, including the shear viscous pressure tensor tt^ u , is written as 



= (e + p)u^u v - pg^ u + vr^ (1) 

where e is the energy density, p is the hydrostatic pressure, and u is the hydrodynamic 4- velocity, 
normalized by u^u^ = 1. T^ v satisfies the energy- momentum conservation law 

8^ v = 0. (2) 

In the Israel-Stewart theory of dissipative fluids jHj, the dissipative fluxes are treated as 
independent thermodynamic variables which satisfy kinetic relaxation equations. The transport 
(relaxation) equations for the shear viscous pressure read 

2 Dtt^ = -—tt^ + V^u^ (3) 

where D = u^d^ is the convective time derivative, rj is the shear viscosity coefficient. is 
related to relaxation time by = 2r//?2- The angular bracket in u v ) is defined as 

V^tt") = -[W + W] - -A^d a u a , (4) 
2 3 

where V M = — u^D is the transverse gradient operator and = g^ v — u^u v is the projector 
orthogonal to the flow velocity u^. 

The viscous pressure tensor tt^ u is symmetric and traceless tt^ v = tt u ^ and tt^ = 0. 
Furthermore, it is transverse to the hydrodynamic 4-velocity, u^tt^ u = 0. Thus it has 5 
independent components. 

There are 10 unknowns (e, p, three components of the hydrodynamic velocity u, and 5 
viscous pressure components) and 9 equations (4 energy-momentum conservation equations and 
5 transport equations for the independent components of tt^ u . The system is closed by the 
equation of state p = p(e). We assume longitudinal boost-invariance and cylindrical symmetry 
in the transverse direction. These additional symmetries reduce the independent variables to 
four (energy density e, radial velocity v r and two components, say ir rr and tt^, of the viscous 
pressure tensor). In (r,r,(f),r)) coordinates, the hydrodynamic and relaxation equations read 
(see appendix) 



d T T TT + d r (y s T^) = - P - l-L , (5) 

d T T^ + d r (v r T Tr ) = P —— d r [p + —y (6) 

d T ir rr + v r d r 7r rr = — fvr rr -27 ? V (r n r)N ) , (7) 

Ttt7 v ' ' 

.— U^-2t]^u^\. (8) 

Ttt7 v 



dq-TT + V r d r TT 4 



where A^ v = rrA^ v , p = rrp and v s = T Tr /T TT '. 

The four equations are solved simultaneously using the SHASTA-FCT algorithm [Tl|. At 
each time step, we obtain the energy density e and radial velocity v r from T TT , T TV , and Tr rr by 
iteratively solving (in analogy to the ideal fluid case |15| ) 

rprr 

e = T TT - v r T Tr . (9) 



T TT + p(e) + 7r rr /7 2 ' 




Figure 1. Contours of constant temperature 
in the r— r plane for ideal and non-ideal 
hydrodynamic evolution, using Tj = 0.3 GeV 
at Ti = 0.5 fm/c. 



Figure 2. Contours of constant temperature 
in the r— r plane for ideal and non-ideal 
hydrodynamic evolution, using Tj=0.5 GeV 
at Ti = 0.1 fm/c. 



3. Equation of state, transport coefficient and initial conditions 

The equation of state is one of the important inputs for hydrodynamics. It is the link between the 
macroscopic and microscopic world. In this exploratory work, we have used a simple equation 
of state, p = \e, with e = aT 4 , a = (16 + ^N f )^. 

In classical kinetic theory, explicit expressions can be obtained for the viscosity coefficient 
rj and the relaxation time r n in terms of the collision term. For a strongly coupled QGP, 
neither r] or r n are known. We treat them as phenomenological parameters. For guidance, we 
use perturbative ^1 EI! an d AdS/CFT ^H] estimates for rj, respectively, and a kinetic theory 
estimate jH] for t^. 

The shear viscosity coefficient rj for hot QCD was determined perturbatively to leading 
logarithmic accuracy in ^JEI- For a s 0.5 the result in ^7] gives 

- = 0.135. (10) 

s 

Recently the shear viscosity was also evaluated in a strongly coupled gauge theory (N = 4 SUSY 
YM theory), using the AdS/CFT correspondence [T^]. The corresponding AdS/CFT estimate 
of shear viscosity is 

- = 0.08 (11) 



In kinetic theory, in the Boltzmann gas approximation, the relaxation time is estimated as 

Tn = 277/32=27? | 0- 

For the initial energy density distribution in the transverse plane, we use the Woods-Saxon 
parameterisation: 

e(r) = —% =w , (12) 
1 + e~ 

with R = QA fm, a = 0.54 fm. This is not very realistic, but facilitates comparison with the 
results of JB]- (eo = oT^ is the central energy density at initial time t = T{.) We have considered 
two sets of initial conditions: (i) Tj = 0.3 GeV and Tj = 0.5 fm/c, and (ii) Tj = 0.5 GeV and 
Tj = 0.1 fm/c. Both sets have similar total initial entropy content. 

We also assume that at the initial time r, the radial velocity vanishes (v r (r) = 0). 

For the non-ideal fluid, initial viscous pressures, 7r rr and ir^ are required. Even though v r 
and its derivatives are zero initially, due to the Bjorken longitudinal motion the stress tensor is 
not zero: V' r it r ' > = r 2 V^u^ > = 3^-. We assume that at initial time Tj, the viscous pressure 
components are fully relaxed to the Bjorken scaling expansion values, 

vr- = rV* = (13) 
3r, ; 




Figure 3. Freeze-out surface 
at T/=160MeV for different 
initial viscous pressures. 



Figure 4. Freeze-out surface 
at T/=160MeV for different 
shear viscosities 77. 



Figure 5. Freeze-out surface 
at T/=160MeV for different 
relaxation times r n . 



4. Results 

4-1- Space-time evolution of ideal and non-ideal QGP 

With have solved the hydrodynamic equation for ideal and non-ideal fluids, with identical initial 
conditions. We do not invoke a phase transition. The initial QGP evolves without undergoing 
phase transition until it freezes out at a an assumed freeze-out temperature Tj = 160 MeV. 

In Figs-^andEl we show contours of constant temperature in the r— t plane. The black lines 
are for the ideal fluid, while the red and blue lines are for non-ideal fluids with 7y/s=0.08 and 
0.135, respectively. For ideal hydrodynamics the almost vanishing slope of the isotherms near 
r = indicates that the transverse rarefaction wave has not reached the fireball center before 
the system decouples at Tj = 0.16 GeV. The cooling at r = is controlled by the longitudinal 
scaling expansion. Consequently, the Bjorken cooling law, T 3 r=constant, is well maintained at 
r = 0. For non-ideal flow, the Bjorken cooling law is violated. Viscosity opposes expansion and 
cooling. Thus the fluid cools more slowly and the lifetime of the fluid is extended. For a fluid 



with Tj = 0.3 GeV at r, = 0.5 fm/c, the lifetime at r = is extended by 30% and 55% for ry/s=0.08 
and 0.135, respectively. For the higher initial temperature (Tj = 0.5 GeV at Tj = 0.1 fm/c), the 
viscosity effects are even more prominent: in the fireball center the fluid life time is enhanced by 
more than 120% and 170% for rj/s=0. 08 and 0.135, respectively. One also observes an increased 
transverse expansion. 

The space-time evolution of non-ideal fluids depends sensitively on the (i) initial viscous 
pressure, (ii) viscosity coefficient and (iii) relaxation time. In Fig. |3] the freeze-out surface at 
TV = 0.16 GeV is shown for different initial viscous pressures, 7r rr (= r 2 7r^)=0, ^- and 4^-, 
respectively, using 77/5=0. 135. The higher the initial viscous pressure, the more extended is the 
freeze-out surface. The life time of the dissipative QGP is extended by 20% if the initial viscous 
pressure is increased from zero to The freeze-out surface also depends sensitively on the 
value of the viscosity coefficient (Fig.Hj). As the viscosity decreases, departure of the freeze-out 
surface from ideal behavior also decreases. In Fig. |S]we show the freeze-out surface for different 
relaxation times, t w = 0.5rk in , Tk in and 2rk in (where 7\ in = t^), for fixed viscosity n/s = 0.135. 
As the relaxation time is increased by a factor 4, the freeze-out time in the fireball center is 
decreased by 25%. 

The evolution of the radial expansion velocity is also changed in viscous fluids. In Fig. [H] 
we show lines of constant radial flow v r in the r— r plane, for the ideal fluid (black line) and 
non-ideal fluids (red and blue lines). The radial velocity grow faster in the non- ideal fluid. 
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Figure 6. Contours of constant radial 
velocity in the r— r plane for ideal and non- 
ideal hydrodynamic evolution, using Tj=0.3 
GeV at Ti=0.5 fm/c. 
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Figure 7. The final parton px distri- 
bution at freeze-out (TV = 160MeV). Initial 
temperature and time are Tj = 0.5GeV and 
Tj = 0.1 fm/c. 



4-2. Viscous corrections to the distribution function 

The changed space-time evolution of the dissipative QGP may have significant effects on particle 
production. Although we cannot compute the pt spectra of final hadrons, since we did not 



include the hadronization phase transition, we can calculate the pt spectra of partons on the 
T = Tf freeze-out surface. 

When viscous corrections are included, the distribution function is modified from its local 
thermal equilibrium form [6*1 IT3]. When the viscous corrections are small, the distribution 
function can be written as 

f(x,p) = f eq (x,p)(l + Af) (14) 
with the equilibrium distribution function 

f ( x v) = — - « JL e -up/T 

If Af is restricted to 2nd order in its form can be determined as pi fTTl IT3] 

Af = C Vil p v ^ v (15) 

where C = 2( J p)T , . 

In Fig. we show the parton px spectra at Tf = 160 MeV. The black line is for the ideal fluid. 
The red and blue lines are for non-ideal fluids with rj/s=0. 08 and 0.135, respectively. Viscous 
flow is seen to enhance high px production by a factor 5-10, depending on viscosity. However, 
since n^ v ^ already at initial time Ti, already the initial px spectra for the viscous fluids are 
flatter than for the ideal fluid. To obtain the same measured final py-spectrum in the viscous 
case thus requires a retuning of initial conditions. This is presently under investigation. 



5. Summary and conclusions 

We have studied the effects of shear viscosity on the space-time evolution of a relativistic 
QGP fluid. In the 2nd order theory of dissipative fluid dynamics, the dissipative fluxes are 
treated as thermodynamic variables which follow kinetic relaxation equations. These are solved 
simultaneously with the energy-momentum conservation equations. Viscosity opposes expansion 
and cooling. Consequently, for a fixed freeze-out temperature, the freeze-out time (fluid lifetime) 
is increased. The evolution of the radial flow velocity is also changed - the viscous fluid generates 
more transverse flow. Viscous effects flatten the pt distribution and require a retuning of the 
initial conditions if the final spectra are to be kept fixed at their measured values. 
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Appendix A. Coordinate transformations 

Instead of Cartesian coordinates = (t, x, y, z) we use curvilinear cylindrical coordinates in 
longitudinal proper time and rapidity, x m = (r, r, (j>, rj): 



t 

x 



t cosh rj; 
r cos (j>; 

r sin <p; 
t sinhr/; 



Vt 2 ^ 2 
x 2 + y 2 



tan 
1 



-ify 



. t + Z 

- In . 

2 t-z 



The differentials 



dt = dr cosh rj + drjT sinhr/, 

dx = dr cos <p — d<pr sin <p, 

dy = dr sin cf> + d(j>r cos <f>, 

dz = dr sinhr] + drjT cosh r/, 



and the metric tensor is easily read off from 



ds 1 



namely 
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In curvilinear coordinates we must replace the partial derivatives with respect to x^ by covariant 
derivatives (denoted by a semicolon) with respect to x m : 



-life 



V I -pi rprak , rpim-pk 

Qj-P P m m P' 



The only non-vanishing Christoffel symbols are 



vv 



TTj 



Tl = l/r; 



-r: 



l/r. 



(A.11) 



The hydrodynamic 4-velocity = ~f(l,v x ,v y ,v z ) is transformed to u m = 7(1, v r , 0,0), with 
7 = v 2 . From here on, we drop the bars over tensor components in x-coordinates for 

simplicity. 



Appendix B. Relaxation equations for the viscous pressure tensor 

Being symmetric and traceless, the viscous pressure tensor ir^ u has 9 independent components. 
The assumption of boost invariance reduces this number by 3 (V^ w u r? '=0, m ^ 
rf). The additional assumption of cylindrical symmetry eliminates two more components 
(V^ m = 0, rj ^ m ^ (p. The transversality condition u m ir mn = eliminates another two 
components (u<p and u v vanish and thus yield no constraints). Thus, with boost-invariance and 
cylindrical symmetry, the viscous pressure tensor has only two independent components which 
we here choose as 7r rr and ir^. 

The viscous pressure tensor relaxes on a short kinetic time scale 77 to 2r/ times the shear 
tensor a^ v = V^u^ The relaxation equations for 7r rr and tt^ are 



d T 7T rT + V r d r 7T Tr = — U TT -27]V {r u r A , (B.l) 

d T ir^ + v r d r 7r^ = —(ir^- 2r}V {,x u®) . (B.2) 



Tn7 

The (rr) and (4>(j)) components of the shear tensor can be written as 



V (r u r) = _Q rU r _ u r Du r _l A rr d . u = 1 2 _Q Tl _ dr ( lv ^\ (R3) 

V<*«*> = V^.u = lf^-^V (B.4) 

r z r 3 r V 3 r J 



with the scalar expansion rate 

d-u = d T u T + 8 r u r + ^ + y = (V + ~) 7 + (d r + ~) (7« r ) (B.5) 
and the convective detivative 

D = u-d = 7(<9 T + v r d r ). (B.6) 

The projectors A rr and A 9 ^ are, 

A rr = 5 rr - u r u r = -1 - (u r ) 2 = - 7 2 

The remaining non- vanishing components of the viscous pressure tensor 7r rT , 7r Tr , and 7r w are 
eliminated using the constraints u m TT mn = and g mn '^ mn = 0, yielding 

1 / TV \ 

7r rr = v r ir Tr = v?ir rr , vr w = --J ^-^- + rV*J . (B.7) 
Appendix C. Energy-momentum conservation 

With longitudinal boost-invariance and cylindrical symmetry, the energy-momentum conserva- 
tion equations T mn - n = yield 



d T f- + dr ( Vs f-) = J-- rr h'-^\ (c.i) 



<V/' rr + M'V' r '') ■ P + rV + ((V2) 



7' 



where i m " = rry4 mn , p = rrp, and u s = T Tr /T TT . 



